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^^ I We present new sets of fragmentation functions for D*^ mesons, both at leading 

and next-to-leading order. They are determined by fitting LEPl data on inclusive 



D*^ production in e+e~ annihilation. In one of the sets, we take the charm-quark 
fragmentation function to be of the form proposed by Peterson et al. and thus 
flj ■ obtain updated values of the ec parameter and the c — > D*^ branching ratio. The 

^^ ■ new fragmentation functions lead to an excellent description of other e^e~ data 

with centre-of-mass energies between 10 and 35 GeV. They also nicely agree with 
recent HERA data on inclusive D photoproduction in ep collisions, which may 
be considered as a test of the universality of the fragmentation into D mesons. 
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1 Introduction 

Recently, the HI ||l| and ZEUS @] collaborations at HERA presented data on the differ- 
ential cross section d?a /dyiah dpr of inclusive D*^ production in low-Q^ collisions, where 
T/iab and pt are the rapidity and transverse momentum of the produced D*^ mesons in the 
laboratory frame, respectively. These measurements extended up to pt = 12 GeV. These 
data were compared with next-to-leading order (NLO) QCD predictions in the so-called 
massless- charm scheme and good agreement was found. 

In the massless-charm scheme, the charm-quark mass rric is neglected, which should be 
a reasonable approximation for pj- ^ rric- In this approach, charm is considered to be one 
of the active flavours in the same way as the lighter u, d, and s quarks. Then, the collinear 
singularities corresponding to the a^ ln{p'^/m'^) terms in a scheme where the charm-quark 
mass is finite and only three active flavours are taken into account, are absorbed into 
the charm-quark parton density functions (PDF's) and into the fragmentation functions 
(FF's) of charm quarks into D*^ mesons. Thus, in a NLO analysis, the following ingredi- 
ents are needed: (i) the hard- scattering cross sections for the direct- and resolved-photon 
processes calculated in the massless approximation {rric = 0) with nj = 4 active flavours 
and with the initial- and final-state collinear singularities, including those of the charm 
quark, subtracted; (ii) the PDF's of the proton and the resolved photon, where charm is 
treated as a light flavour; [iii) the FF's characterizing the hadronization of the massless 
partons, including the charm quark, into D*^ mesons. (Similarly, the bottom quark is 
treated as a massless parton above a certain threshold, leading to n^ = 5.) 

This massless approach was originally considered in [Q as a possibility to make pre- 
dictions for large-pr heavy-quark production, and it was first applied to the production 
of large-pT hadrons containing bottom quarks in pp collisions 0. Subsequently, it was 
employed to study charm-quark production in ep 0, 7p 0, and 77 collisions p. 

In order to obtain definite predictions for D*^ production cross sections in the massless- 
charm scheme, one obviously needs the PDF's of the charm quarks in the proton and the 
photon, which do not enter the massive calculation, as well as the FF's for D*^ mesons. 
Both the PDF's and FF's are basically non-perturbative input and must be determined 
by experiment. 

This offers the possibility to use the data of D*^-photoproduction experiments in re- 
gions where the massless-charm approach is a good approximation, to obtain information 
on the charm PDF's of the proton and the photon and on the FF's of the quarks, in 
particular the charm quark, and the gluon into D*^ mesons. 

The FF's of the D*^ mesons can be determined from measured cross sections of D*^ 
production in e~^e~ annihilation. In these experiments, the centre-of-mass (CM) energy 
^/s of the electron-positron system is large compared to rric, so that the massless-charm 
approach is well justified in this case. Then, one is left with the task of constraining the 
charm components of the proton and the photon. The photon PDF's enter, for instance, 
the resolved-photon cross section of inclusive D*^ photoproduction in ep collisions. In 
a previous work by two of us with Spira 0, it was found that, at moderate pt, the 
resolved-photon cross section is of the same order of magnitude as the direct-photon cross 



section. Furthermore, the additional charm component in the proton turned out to only 
contribute marginally to the cross section of inclusive D*^ photoproduction. Therefore, 
the study of this cross section offers an ideal possibility to specifically learn about the 
charm component of the photon, on which not very much is known from other experiments. 
The traditional place to investigate the charm PDF of the photon has been fully inclusive 
deep inelastic charmed-hadron production in 67 collisions. Unfortunately, very little is 
known experimentally so far [p. 

Several NLO parameterizations of the charm component of the photon exist in the 
literature [llO| , |rTi [T2| , and we shall use them for our predictions in the following. The point 



to be made is that these charm PDF's have not yet really been tested against experimental 
data which are particularly sensitive to it. 

In the previous work on D*^ production |^, we considered several approximations for 
the fragmentation of charm quarks into D*^ mesons, including also the Q^ evolution to 
higher scales. The most realistic description was by the Peterson FF |]T3| with evolution. 



which was also used for the comparison with the HI [Q] and ZEUS [g] data. Apart from 
the overall normalization, the distribution of Peterson et al., which is used as input at 
the starting scale /iq, depends only on one parameter, ec- The value of e^. was taken from 



the work by Chrin |]I4[, which is a phenomenological analysis of data on charmed-meson 
production (not necessarily restricted to D*^ production) in the PETRA-PEP energy 
range of e^e~ annihilation. In the meantime, much better data, specifically on D*^ 
production, have been made available by the ALEPH ||T5[ and OPAL |T^ collaborations 



at LEPl. These data should give us an excellent handle on the Z)*^ FF's, so that we need 
not rely on Chrin's value for the ec parameter any more. 

It is the purpose of this work to improve the description of the charm and bottom FF's 
into Z)*^ mesons by fitting to the ALEPH and OPAL data on the basis of NLO evolution 
equations. These FF's will be tested against D*"^ production data at lower CM energies 
13 coming from PETRA JTH, PEP JTl, and DORIS |20|. They will also be used to make 



theoretical predictions for inclusive D*^ photoproduction at HERA, following the earlier 
work 0. 

The outline of this paper is as follows. In Section 2, we shall briefly describe the NLO 
formalism of charm-quark fragmentation including the transition from massless to massive 
factorization of final-state coUinear singularities also employed in [Q. In this section, we 
shall specify the D*^ FF's at the starting scale /xq, adopting two different functional 
forms. We shall consider not only charm-quark fragmentation, but also bottom-quark 
fragmentation which is needed to describe the LEPl data. Finally, we shall present the 
resulting NLO predictions for D*"^ production in e^e~ annihilation and compare them 
with data from ALEPH and OPAL used in the fits, and with the other e^e~ data at lower 
CM energies. In Section 3, we shall confront our predictions for ep -^ D*"^ + X based 
on the two forms of FF's with the ZEUS data 0. Our conclusions are summarized in 
Section 4. 



2 D*^ Production in e^e Annihilation 

To construct LO and NLO sets of FF's for D*^ mesons, we make use of the new LEPl 
data from ALEPH ]15|] and OPAL ]16|. In hadronic Z-boson decays, charmed mesons 



are expected to be produced either directly through the hadronization of charm quarks 
in the process Z ^ cc or via the weak decays of B hadrons produced in Z — > bb, with 
an approximately equal rate. Charmed mesons from Z ^ cc allow us to measure the 
charm-quark FF. The main task of both experimental analyses was to disentangle these 
two main sources of D*^ production at the Z-boson resonance. In the two papers p!5|JT6| 



a number of different techniques are used to tag D*^ mesons from B decay. The D*^ 
mesons are reconstructed in the decay chain D*^ -^ D^n^ and D^ —>■ K~tx^ , which allows 
for a particularly clean signal reconstruction. 

The results of the two measurements are separate, differential yields for D*"^ mesons 
from Z ^f cc and Z ^ bb decays, normalized to the total number of multihadronic Z- 
boson decays, as a function of x = 2E{D*^)/ y/s, where E{D*^) is the measured energy 
of the D*^ candidate. The contribution due to charm-quark fragmentation is peaked at 
large x, whereas the one due to bottom-quark fragmentation has its maximum at small 
X. For the fitting procedure, we use the x bins contained within the interval [0.1,0.9] 
and integrate the theoretical functions over the bin widths, which is equivalent to the 
experimental binning procedure. As in the experimental analyses, we sum over D*^ and 
D*^ mesons. As a consequence, there is no difference between the FF's of a given quark 
and its antiquark. 

We take the starting scales for the D*^ FF's of the gluon and the u, d, s, and c quarks 
and antiquarks to be /xq = 2mc, while we take /iq = 2mh in the case of bottom. The 
FF's of the gluon and the first three flavours are assumed to be zero at the starting scale. 
These FF's are generated through the /i^ evolution, and the FF's of the first three quarks 
and antiquarks coincide with each other at all scales /i. For the parameterization of the 
charm- and bottom-quark FF's at their respective starting scales, we employ two different 
forms: 

(i) In our standard set (S), the FF's of the charm and bottom quarks are both assumed 
to be of a form usually used for the FF's of light hadrons, namely, 

DQix,fxl)=Nx''{l-xf, (1) 

where Q = c,b. Here and in the following, we skip the index D*^, since we only consider 
D*^ production throughout this paper, e.g., Dc denotes the FF of the charm quark into 
D*^ mesons, etc. 

(ii) In the so-called mixed set (M), we use for the bottom-quark FF the form of (|I|) 
and for the charm-quark FF the Peterson distribution [J13i, 



The Peterson form is particularly suitable to describe FF's that peak at large x. Since 
the bottom-quark FF has its maximum at smaller x values, we obtain intolerably bad fits 



if we also use (^ in the case of bottom. This is the reason for choosing the bottom-quark 
FF at fiQ = 2mh in set M to be of form (|l|). 

Obviously, the functional form (|l]) has more flexibility, since, apart from the normal- 
ization factor N, it depends on two parameters, a and (3, whereas there is just a single 
parameter, ec, in the Peterson distribution. 

The cross section of inclusive D*^ production in e~^e~ annihilation. 



e+e- -^ (7, Z) -^ D*^ + X 



(3) 



is calculated from the same formula as for any other hadron, 
1 d(r{e+e- ^ D*^ + X) 



C^tot 



dx 



O"tot a J^ ^ V-2 ■> J dz ^ ■' ^ 



(4) 



Here, the sum extends over all active partons {a = g,u,u, ... ,b,b), ^ is the renormalization 
scale of the partonic subprocesses, and Mf is the fragmentation scale. At NLO, Mf defines 
the point where the divergence associated with coUinear radiation off parton a is to be 
subtracted. This applies to the u, d, and s quarks as well as for the c and b quarks, 
the masses of which are neglected. To NLO in the modified minimal subtraction (MS) 
scheme, the cross sections of the relevant subprocesses are given by 
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Here, Cg. are the effective electroweak couplings of the quarks g^ to the Z-boson and 
the photon including propagator adjustments, and P^^i, are the LO timehke splitting 
functions, which may be found in the Appendix. The C functions are given by [^ 
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[21nx + ln(l ~x)], 



(6) 



where C(2) = 7r^/6, Cp = 4/3, and the plus distributions are defined in the Appendix. We 
evaluate a^ using the two-loop formula with nf = 5 flavours. We identify fi = Mf = y/s, 
so that, in (|^), the terms involving ln(s/M|) vanish. 

Up to this point, the formalism is identical to the one for describing the fragmentation 
into light hadrons based on the MS factorization scheme. For the fragmentation into heavy 
particles, we take one further step. As in the previous work [^, we adjust the factorization 
of the final-state coUinear singularities associated with the charm and bottom quarks in 



such a way that it matches the massive calculation. We achieve this by changing the 
factorization scheme. Specifically, we substitute in the hard-scattering cross sections (||) 
o^s{^^')P:i'i\x) Hs/Mj) -. a,{f^')P^STi^) Hs/Mj) - a, (/i^) dQ,{x), where 



dQQ{x) = -Pr 



y(0,T), 



/^O 



x)\n^ + CF<-26il-x)+2 
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X 



+ 4 



ln(l — x) 
1 — X 



- (l + x)[l + 21n(l-x)]} 
dgqix) = dQg{x) = dQQ{x) = 0, 



The (ig^ functions emerge by charge conjugation 
p2| . Note that dQg{x) and dQg{x) do not yet enter 



(7) 



These 



with Q = c,b and q = u,d, s. 
functions are extracted from 
at NLO. The same substitutions must also be performed in any massless hard-scattering 
cross sections that are to be convoluted with our FF's, in particular with those of inclusive 
particle production in ep collisions, which will be studied in the next section. This change 
of factorization must be regarded in connection with our assumptions concerning the FF's 
at the input scale, which refer to massive charm and bottom quarks. This means that, 
without this change, we would proceed as in the case of fragmentation into light hadrons 
and obtain somewhat different parameters for the input distributions. 

The FF's in (^) are parameterized at the starting scale Mf = /iq by the ansatze (^) 
and/or (|l]), depending on whether set M or S is considered. Once we know the FF's at 
the scale /iq, the /x^ evolution is ruled by the Altarelli-Parisi (AP) equations |2^. Our 
task is thus to determine the parameters of the distributions (^ and (0) at the starting 
scale /io in such a way that the evolved FF's, when convoluted with the hard-scattering 
cross sections, fit the ALEPH |T^ and OPAL \W\ data at Mf = Mz- In the case of 
form dll), the AP equations can be easily solved with the help of the Mellin-transform 
technique, since the moments of the starting distribution (|1]) are simple. Further details 
are given, for example, in our earlier work [Q. In the case of the Peterson form (|]), the 
computation of the moments in Mellin space is numerically rather complicated; see Q] 
for the corresponding formulae. In this case, we find it more convenient to perform the 
evolution of the FF's in x space by iteratively solving the AP equations in their integral 
form, 



D,{x,t) = D,{x,0) + I dt'Y.1 -PS ( -,a,(t') ) D,{z,f) 

Jo ^ Jx Z \Z J 



(8) 



where t = In [fi"^ / fi^) . The first term on the right-hand side of (§) is the distribution 
at the starting scale /iq. As a check for the numerical stability of our procedure, we 
compared both methods of evolution at LO and NLO for forms (|T]) and (^, and found 



good agreement. The NLO timelike splitting functions P^ 



(i,T), 



X) 



which enter (|^), may be 



extracted from the literature 25 



Since, to our knowledge, there is no reference where all 
these functions are presented in a convenient form, i.e., with the coefficients of the delta 
functions and plus distributions explicitly displayed, we collect them in the Appendix. 
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In the case of set S, we need to fit A^, a, and (3 in (|I]) for the charm- and bottom- 
quark FF's, while, in the case of set M, we need to determine these parameters for the 
bottom-quark FF as well as A^ and ec in (|^) for the charm-quark FF. Both ALEPH [15] 
and OPAL ||l6l extracted the total D*^ momentum distribution as well as the individual 
contributions due to the Z ^ cc and Z —>■ bb channels. In the case of ALEPH, we received 



the individual cc and bb data sets in numerical form via private communication |26|. In 
the ALEPH analysis, the branching ratio of D^ -^ K~tt'^ (or the charge-conjugated final 
state) was taken to be 0.0371. Since this value has been superseded by the value 0.041 p7|. 



which was also used in the OPAL analysis, we rescale the ALEPH data with the factor 
0.9049. In our LO (NLO) analysis, we use the one-loop (two-loop) formula for as with 
A^ = 108 MeV (227 MeV) as in our previous work on inclusive pion and kaon production 



24 1 . At first sight, the value A^ = 108 MeV of our LO fit may seem unrealistically low; 



MS 



this is, however, inconsequential because a higher value can easily be accommodated by 
slightly modifying the other fit parameters without damaging the quality of the fit. The 
starting scales are chosen to be /iq = 2mc with rric = 1.5 GeV for all FF's except for the 
case of bottom, where it is taken to be /iq = 2mf, with m^ = 5 GeV. 

The values of A^, a, (3, and ec resulting from the combined LO and NLO fits to these 
data are displayed in Table |I|. The parameter ec appears only in the case of the charm- 
quark FF's of set M. In the case of the charm-quark FF's of set S, we have a ^ (3, 
as is expected for a distribution peaked at large x. On the other hand, /5 > a for the 
bottom-quark FF's, which have their maxima at smaller values of x. The parameters of 
the bottom-quark FF's of sets S and M are very similar, i.e., the data fix the charm- and 
bottom-quark FF's independently of each other. 

In Figs. |l|a and b, we compare the ALEPH data [p!5|j26|] with our fits according to sets 



S and M, respectively. In Figs. ^ and b, we do the same for the OPAL data [Q. In each 
case, we show both LO and NLO results. Except for very small x, the LO and NLO curves 
are very similar. At small x, we observe large differences between LO and NLO, indicating 
that, in this region, the perturbative treatment ceases to be valid. In this region, also the 
massless approximation is not valid anymore. This may be understood by regarding the 
phase-space boundary for the production of D*^ mesons. Since D*^ mesons have mass 
m{D*^) = 2.11 GeV [|^, they can only be produced for x > Xmin = 2m{D*^)/y/s. At 
a/s = Mz, one has a;min = 0.046. This is just the region where our NLO results turn 
negative. Therefore, our results are meaningful only for x^Xcut = 0.1. Consequently, 
we exclude the first bin of the OPAL data (0 < x < 0.1) from our analysis. Note that 
this phase-space boundary is not only theoretically motivated, but does also turn up in 
the measurements, which give values consistent with zero whenever x < Xmin- This is not 
only true for the OPAL data, but also for data taken at lower energies, e.g., those from 
ARGUS O. 

In Figs. |]a, b and ||a, b, we have plotted the fully normalized cross section of (^), which 
also includes the gluon contributions in NLO. The shown charm and bottom contributions 
include terms with a = c,c and a = 6,6 in (^, respectively, as well as the term with 
a = g distributed to the charm and bottom terms in proportion to e"^./ Y^^lie"^,, which 



is motivated by (^. The contributions from a = u,u, d, d, s, s are negligible in size; they 
add up to less than 2% of the total integrated cross section and mostly contribute at small 
X. This procedure approximately produces the quantities that should be compared with 
the ALEPH and OPAL data. For the time being, we ignore the effect of electromagnetic 
initial-state radiation, which has not been corrected for in the data. 

The quality of the fit is measured in terms of the average x^ for ^ill selected data 
points. In Table H, we list the x^ values obtained with sets S and M at LO and NLO 
for the ALEPH and OPAL data, and their combination. We also specify the individual 
results for the Z ^ cc and Z ^ bb samples of ALEPH and OPAL. We observe that the 
Z ^ cc samples yield worse x^ values than the Z ^ bb samples. The OPAL data tend 
to give better x^ values than the ALEPH data. Especially in the case of ALEPH, set S 
leads to better x^ values than set M, which is not surprising if one recalls that set S has 
six degrees of freedom while set M has only five. 

The integrals of the charm- and bottom-quark FF's into D*^ mesons over x give the 
branching fractions of these transitions. For the reasons given above, we restrict our 
considerations to the region Xcut < 2; < 1, with Xcut = 0.1, and define 

5q(^)= /' dxDQix,fx^), (9) 

where Q = c,b. Experimentally, the contribution from the omitted region < a; < Xcut is 
close to zero with a large error. It is interesting to study Bq^jj) for Q = c,b aX threshold 
/i = 2mQ and at the Z-boson resonance /i = Mz- In Table ^, we present the results for 
sets S and M at LO and NLO. We observe that these branching fractions are changed 
very little through the evolution from 2mQ to Mz-, and they are almost equal for charm- 
and bottom-quark fragmentation. The branching fractions at the Z-boson resonance can 
be compared with the experimental numbers reported by the OPAL collaboration |[T6| , 
which are Bc{Mz) = 0.259 and Bb{Mz) = 0.252, with errors of approximately 25%. These 
results, which are obtained from independent fits to all experimental data, are consistent 
with our fits. At LO (NLO), the charm to bottom ratio of the branching fractions at /x = 
Mz is 0.92 (0.94) for set S and 1.07 (1.01) for set M. This agrees well with the experimental 
results B^{Mz)/BbiMz) = 1.15+J-J°(stat.) |T§ and 1.03 ± 0.11(stat.) ± 0.10(syst.) |TB|. 



Another quantity of interest is the mean momentum fraction, 

1 /"i 
(a;)Q(/i) = p . ^ / dxxDqix^iJ,^), (10) 

where Q = c,b. In Table H, we collect the values of {x)q{^) for Q = c,b evaluated with 
sets S and M in LO and NLO at /i = 2mQ, Mz- The differences between sets S and M and 
between LO and NLO are minor. However, the effect of the /z^ evolution is significant, 
e.g., {x)c{fJ') moves from approximately 0.7 at /i = 2mc down to about 0.5 at /x = Mz- 
Our results for {x)c{Mz) should be compared with the experimental numbers reported by 



ALEPH fl5| and OPAL M 



{x),{Mz) = 0.495 ^°'° J° (stat.) ±0.007(syst. 



(a;),(Mz) = 0.515 ^°'°°^stat.) ± 0.010(systO 
— U.UU5 



:iii 



respectively. We conclude that our fits are in reasonable agreement with the independent 



determinations in ||T5|,|16|. 

Next, we test whether the FF's constructed from the ALEPH and OPAL data, after 
evolution to lower scales, lead to predictions for D*"^ production by e^e~ annihilation that 



are consistent with experimental data taken at lower CM energies, 
surements at lower energies were performed, e.g., at DORIS, CESR, 
for a compilation, see 



A number of mea- 
PEP, and PETRA; 
These data have much lower statistics than the data from 
LEPl. For the comparison, we selected the data from ARGUS [^ at ^/s = 10.49 GeV, 

at -y/i 



m- 



from HRS 19 



at ^/s = 34.2 GeV. 
meson, namely, D^ — 



The TASSO 
K- 



TT TT TT 



29 GeV, and from TASSO 
collaboration measured two decay channels of the D^ 
and D^ —>■ K^n^. The latter decay channel also allowed measurements at lower values of 
x; we shall denote this data set by TASSO 2 and the other one by TASSO 1. For consis- 
tency, we corrected the HRS and TASSO data by updating the values of the branching 
ratios of D*^ -^ D°7r+, D° -> K-n+, and D^ -^ R-n+Tr+n- according to the 1994 
tables of the Particle Data Group p^, which were also used by OPAL [|1^. This leads 
to correction factors of 0.702 for TASSO 1 and 0.677 for TASSO 2. In the publication 



by HRS [0, the cross sections are not given separately for the different decay channels, 
so that we can only approximate the required correction by an overall factor of 0.877. 
In Fig. ^, we compare the ARGUS, HRS and TASSO data on the cross section sda/dx 
of e+e" — *> D*^ + X with our respective LO and NLO predictions based on set S. The 
theoretical results are calculated according to (^) with nj = 5 quark flavours (except 
for the case of ARGUS, where nj = 4 is used), i.e., also the contributions due to gluon 
and light-flavour fragmentation are consistently included. We only show our results for 
0.1 < X < 0.9. For completeness. Fig. ^ also includes the comparison of the OPAL data 
T6| on the dimensionless cross section {l/atot)da/dx with our LO and NLO predictions 



according to set S. As can be seen from Fig. |^, the agreement between our predictions and 
the data is satisfactory. This nicely demonstrates that the data indeed exhibit the scaling 
violation predicted on the basis of the LO and NLO AP evolution equations for the FF's. 
In fact, the change in the shape of the differential cross section with /i^ is mainly due to 
the bottom quark. The ARGUS data |^ are taken off the resonance, at y/s = 10.49 GeV, 
where the bottom quark is not yet active. The inclusion of the bottom-quark contribution 
leads to a softening of the distribution, as may be seen in the case of the HRS data [|TU| at 
y/s = 29 GeV and the TASSO data [|18l at ^i = 34.2 GeV. The evolution from 29 GeV to 
34.2 GeV has no discernible effect. Going to the Z-boson pole at y/s = 91.2 GeV further 
increases the relative importance of the bottom-quark fragmentation and leads to an even 
softer spectrum. For 0.2 < x < 0.8, there is little difference between the LO and NLO 
cross sections. 

The successful comparisons in Fig. |] make us confident that our D*^ FF's, although 
constructed at y/s = Mz, also lead to useful descriptions of D*"^ fragmentation at other 
scales. In the next section, we shall exploit this property together with the universality 
of fragmentation to make predictions for inclusive D*^ photoproduction at HERA. 



9 



3 D*^ Production in Low-Q"^ ep Collisions 

In this section, we compare our NLO predictions for the cross section of inclusive D*^ 
photoproduction in ep scattering at HERA with data from the ZEUS collaboration |0. 
The calculation of this cross section proceeds as in the previous work 0, except that we 
now use the FF's constructed in the last section. For the calculation of the cross section 
(i^cr/t/yiab dpT, we adopt the present HERA conditions, where Ep = 820 GeV protons 
collide with E^ = 27.5 GeV positrons in the laboratory frame. We take the rapidity to be 
positive in the proton flight direction. The quasi-real photon spectrum is described in the 
Weizsacker- Williams approximation by the formula given in 0. This spectrum depends on 
the photon-energy fraction, x = E^/E^, and the maximum photon-virtuality, Q^ax- I^ ^^^ 
ZEUS experiment [0], where the final-state electron is not detected, Q^ax — ^ GeV^ and 
0.147 < a; < 0.869, which corresponds to jp CM energies of 115 GeV <W< 280 GeV. We 
adopt these kinematic conditions in our analysis. We work at NLO in the MS scheme with 



Uf = 4: flavours. For the proton and photon PDF's we use set CTEQ4M |^ with A^ 



296 MeV and set GRV HO [^ converted to the MS factorization scheme, respectively. 

L^l^ TheAill 

'■MS' ^ MS 



We evaluate Q;s(yU^) from the two-loop formula with this value of Ai-=^. The Ai^ values 



implemented in the D FF's and photon PDF's are 352 MeV, which corresponds to 
the value A-^ = 227 MeV quoted in the previous section, and 200 MeV, respectively. 
We identify the factorization scales associated with the proton, photon, and D*^ meson, 
and collectively denote them by Mf. We choose the renormalization and factorization 



scales to be /i = m^ and Mj = 2mT, respectively, where m^ = Jpr + rnl is the D*^ 
transverse mass. Whenever we present LO results, these are consistently computed using 
set CTEQ4L [g|] of proton PDF's, set GRV LO |g of photon PDF's, the LO versions 



of our sets S and M of D*^ FF's, the one-loop formula for a^ with A|-=^ = 236 MeV, and 
the LO hard-scattering cross section. 

The photoproduction cross section is a superposition of the direct- and resolved-photon 
contributions. In our NLO analysis, the resolved-photon contribution is larger than the 
direct one for moderate pt- This statement depends, however, on the factorization scheme; 
only the sum of both contributions is a physical observable and can be compared with 
experimental data. The bulk of the resolved-photon cross section is due to the charm 
component of the photon PDF [Q . 

We first consider the pt distribution da/dpT integrated over the rapidity interval 
—1.5 < yiab < 1 as in the ZEUS analysis 0. In Fig. ^, we compare the respective 
ZEUS data with our LO and NLO predictions based on sets S and M. The theoretical 
results are multiplied with the factor 1.06 to account for the fact that OPAL |]TB| used 



the D*^ and D° branching fractions of [^, while ZEUS used those of |2^. The 



NLO distributions fall off slightly less strongly with increasing px than the LO ones. It is 
remarkable that there is hardly any difference between the set-S and set-M results; in fact, 
the LO-S and LO-M (NLO-S and NLO-M) curves almost coincide. This means that the 
details of charm-quark fragmentation are very tightly constrained by the LEPl data, and 
that the considered variation in the functional form of the charm-quark FF at the starting 
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scale has very little influence on the pt distribution. The NLO curves should be compared 
with Fig. 8c in , where the charm-quark FF at the starting scale /io = 2m(. = 3 GeV was 
taken to be of the Peterson form with ec = 0.06 and Bc{2mc) = 0.260. This corresponds 
to A^ = 0.163. We recall that, in the case of set NLO M, we have e^ = 0.0204, A^ = 0.0677, 
and Bc{2mc) = 0.232; see Tables |I] and ||. We observe that the result for /x = M//2 = m^ 
in Fig. 8c of is somewhat smaller than the NLO results of Fig. ^, but has a similar 
shape. The agreement with the data is good, even at small pt, where the massless 
approach ceases to be valid. 

In Fig. ^, we compare the yi^h line shape of dP a / dyi^h dpT , integrated over the Pt 
interval 3 GeV < pr < 12 GeV, with the respective ZEUS data points 0. As in Fig. |a, 
we consider the four cases LO S, NLO S, LO M, and NLO M. Again, the LO-S (NLO- 
S) results almost perfectly coincide with the LO-M (NLO-M) results, i.e., the functional 
form of the charm-quark FF at the starting scale is of minor signiflcance. The agreement 
between theory and experiment is reasonable. If we compare the NLO results of Fig. ^ 
with the result for n = Mf/2 = m^ in Fig. 9c of 0, we see that the latter is somewhat 
smaller, but has essentially the same shape. 

In [0, the cross section d'^a /dyiat, dpr was also measured at three different values of W. 
To this end, yy^h and px were sampled in the intervals —1.5 < yiab < 1 and 3 GeV < pt < 
12 GeV, respectively. In Fig. ^, we show the W dependence of d'^a /dyigjo dpr, integrated 
over these yiah and pt intervals, and compare it with the corresponding ZEUS data p. 
Similarly to Figs, ^a and b, we consider sets S and M at LO and NLO. Within the rather 
large experimental errors, all four predictions agree reasonably well with the data. In 0] , 
such a flgure was not presented, but it can be found in the ZEUS publication 0]. Again, 
we observe that the input assumptions, encoded in sets S and M, have practically no 
influence on the result. 

When we fltted to the LEPl data, we also had to take into account D*^ production 
by bottom-quark fragmentation, as it gives a sizeable contribution at the Z-boson pole. 
In ep collisions, however, the bottom quark does not contribute at all to inclusive D*^ 
photoproduction below the threshold at mx = mf,. In contrast to e~^e~ annihilation, the 
bottom quark does not contribute signiflcantly even above the threshold. To elucidate 
this point, we repeat the analysis of Fig. |a for Uf = 5 with and without the bottom-quark 
FF and plot the ratio in Fig. ^. We see that, in all four cases, the relative contributions 
due to the bottom-quark FF's are below 5% and decrease with increasing pr- This may 
be understood by observing that, in the case of D*"^ photoproduction in ep collisions, 
the main contribution to the cross section arises at large x, typically at x^O.5, where 
Db{x, M|) ^ Dc{x, M|) because the bottom-quark FF's have a soft spectrum. For the 
same reason, also the contributions of the u, d, and s quarks are negligible. 

We have seen that our theoretical predictions for D*^ photoproduction at HERA are 
very insensitive to the speciflc form of the input distribution of the charm-quark FF, and 
that the influence of the other quark flavours is marginal. When the FF templates of 
sets S or M are fltted to the most accurate data, we obtain an almost unique prediction 
for the charm-quark FF. In fact, the largest remaining uncertainty in the calculation of 
the D*^ photoproduction cross section is related to the charm content of the resolved 
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photon. This offers us the possibihty to employ our well constrained FF's to extract 
information about the charm distribution in the photon. It is well known that the most 
sensitive distributions to study such dependences are those differential in rapidity. This 
fact has been previously exploited, e.g., in |^, where we emphasized the potential to 



determine the gluon PDF of the photon by accurately measuring the inclusive charged- 
particle rapidity spectrum in photoproduction. To illustrate the sensitivity to the charm 
PDF of the photon, we repeat the NLO-S analysis of Fig. ^ with this PDF switched off, 
and compare the outcome with the full calculation in Fig. ^. For comparison, also the 
direct-photon contribution is shown. We observe that 73% to 100% (50% to 89%) of the 
resolved-photon (full) cross section is induced by the charm content of the photon, which 
is most important in the backward direction. 

Inspired by this observation, we now repeat the NLO-S analysis of Fig. ^ using in 
turn the NLO photon PDF sets of 0Jl|, which we denote by GS HO and ACFGP, 



respectively. The authors of |T2[ considered both massless and massive charm quarks. 
In their massive-charm set, the charm-quark PDF is implemented with the boundary 
condition of [^ at threshold. We observe a significant variation of the cross section. 



especially in the backward region, where our predictions for set GRV HO ||T0| and ACFGP 
(massive c) [Q differ by more than a factor of two! Unfortunately, the data are not yet 
accurate enough to rule out any of the available sets. Of course, in order to reliably 
constrain the charm PDF of the photon, it would be preferable to have data with an 
increased low-pj- cut, since our massless charm-quark theory is more reliable in the large- 
Pt region. However, this is a very interesting measurement that has a vast potential for 
improving our knowledge about the charm PDF of the photon. 

4 Conclusions 

In this paper, we determined new sets of D*^ FF's at LO and NLO from fits to precise 
cross-section data obtained in e~^e~ annihilation by ALEPH and OPAL at LEPl. Our 
analysis was based on the massless- charm factorization approach recently developed in |^ . 
In these LEPl data, the fragmentation of charm and bottom quarks into D*^ mesons was 
disentangled, so that the respective FF's are separately constrained. For the description 
at the charm-quark FF, we used two different ansatze at the starting scale, namely, the 
three-parameter standard form conventionally used for light-meson production and the 
two-parameter Peterson form. It turned out that both prescriptions lead to excellent 
fits of the LEPl data, and that the charm- and bottom-quark FF's are each tightly 
constrained. With this method, the fragmentation at any other scale is predictable via 
the AP evolution, which we performed both at LO and NLO. This enabled us to check our 
FF's against e~^e~ data at lower CM energies. In particular, we found very good agreement 
with data from TASSO [18 , HRS [|19|, and ARGUS |20|, which demonstrates that the 



scaling violation encoded in the AP equations is consistent with the data. We stress 
that the results for the starting parameters are correlated with the starting scale. In the 
case of the Peterson description for charm-quark fragmentation at LO (NLO), the epsilon 



12 



parameter and c -^ D*^ branching fraction at the starting scale /zq = 2?7ic were found 
to be ec = 0.0856 (0.0204) and 5e(2m,) = 0.246 (0.232), respectively. This corresponds 
to an X mean value of {x)c{2mc) = 0.644 (0.754). While Bc{ix) is relatively stable under 
evolution to yU = M^, {x)c{lj) is appreciably lowered, down to {x)c{Mz) = 0.490 (0.556). 
We emphasize that the D*"^ FF's thus obtained are universal and may be used to 
make predictions for inclusive D*^ production in any other high-energy experiment. We 
exploited this fact by calculating the cross section of inclusive D*"^ photoproduction at 
HERA, which we compared with the latest data from the ZEUS collaboration. We found 
that the bottom-quark contribution is marginal, and that the results do not depend 
on whether we use the standard or Peterson-type charm-quark FF's, which shows that 
we have good control over the fragmentation aspect of D*^ photoproduction. Under 
the condition that more precise experimental data, in particular in the large-p^ region, 
will be collected at HERA, this will allow one to constrain other non-perturbative input 
information such as the charm-quark PDF of the photon, which is still poorly known. 

NOTE ADDED IN PROOFS 



After the completion of our manuscript, a related paper ||3^ has appeared. The authors 



of p2[ determine D*^ FF's by separately fitting e~^e data from DORIS and LEPl. In 



contrast to our procedure, they also include the dQa functions of (|^ in the evolution, and 
treat the bottom quark on the same footing as the u, d, and s quarks, i.e., they only take 
the charm-quark FF to be non-vanishing at the starting scale. 
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A Appendix: NLO timelike splitting functions 

Although the NLO timelike splitting functions have been calculated by various groups 



25| , we are not aware of any paper where they are presented in a way ready to use, i.e., 
with the coefficients of 6{1 — x) and 1/(1 — x)_|_ displayed explicitly. For the reader's 
convenience, we shall do this here. 

Let Dg{x,jj?), Dq.^x^jj?), and Dq^{x,^'^) be the FF's of the gluon g and Uf quarks g, 
and antiquarks (ji {i = 1, . . . , nf), respectively, into some hadron with momentum fraction 
X at fragmentation scale /x. The /i^ evolution of these FF's is conveniently formulated for 
the combinations 

dI^^ (x, ^^) = Dg^ (x, ^^) ± Dq^ {x, /i^) , 
D^^\x,f,')=Y.Dl^\^,f^') (12) 
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as 
dii^ 



(^A^+)(x,/.2)--iDW(x,/i2)j 



^ ^^ td{T) (X t 2 



(13) 






•J X /C/ \ '^ / 



d^j? 
jj'^d 
diJ? 
fi'^d 
diJ? 



"1 (i^; 



I^^^^a:,/.^) = /^ ^ [Pil), (?"^(^')) ^^^H^,/^^) + 2n,p(!:), (?"^(^'0 ^^(^'^ 






2; 



^pS(p«.(^=))fl"(.,,.^) + pi!;!,(p 



as(/i^)j Dg{z,^^) 



where P^^^fx, Q;s(/i^) j are the timehke a ^ h sphtting functions. PL) (ic, ^^(/i^) j and 

Pg^qfa^, tts(/^^)) are decomposed into valence (non-singlet) and sea (pure-singlet) compo- 
nents as 

i^^J(x,a,(/i2)) =Pj^iP(x,a,(/i2)) ±P,^iP(x,a,(/i2)), 

n^^U^,«.(/"')) =<?(^,«s(^')) +2n^P/iV(^,«.(/^'))- (14) 

All these splitting function have perturbative expansions of the form 

pS, (-.".(^^)) = ^i-Sflx) + (^P^ P^'i^+Oioi) . (15) 



The LO coefficients read 



q-^q 



!„l_,) + 2(_l_)^_l_. 



p;T'w = «'''w = o, 



q-^q 
PitVi^)=CFPFG{x), 

P^tVi^) = ^TnfPGF{x), 



p^tVi^) 



^CA-lTnf]5{l-x)+2CA 



1 — X 



-2 + x-x^ 



X 



(16) 



and the NLO coefficients read 
"3 



p:i7\-)=\cf 



-3C(2) + 6C(3) 



+ CfCa 
U(1-x) + |CfCa 
-|- Cp < —5(1 — x) — -(7 + 3a:) \nx + -(1 -|- x) \n^ x + pFFix) \nx 



M«2' 



'ij + yC(2) - 3C(3)' 


- CpTuf 


[?-2C(2,] 




■-} (i^). 
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--21na; + 21n(l-a;) 



+ CpCa {y^(53 - 187a;) + (1 + a;)[C(2) + Inx] 



+ 7^Pff{x) Ina; f — + Ina; ) |> + CpTuf 



2 2 

-(-1 + llx) - -Pff{x) Inx 

y o 



PqSY\x) = Cf {2Cf - Ca) [2(1 -x) + {l + x)\Yix + Pff{-x)S2{x)] 



p,'iY\-) 



CfT 



2 + a; H x^] - {b + 9x+ -x^] \nx + (1 + x)\n^ x 

9a; 9 / V 3 



P^P(a;) = C^ |1(-1 + 9x) + (^-8 + I) Ina; + 2a;ln(l 



-x)+[l 



X 



\v? x + pfg{x) 



-8C(2) + 41nxln(l - x) + \n^{l - x) - 85i(a;)]} + CfCa |^ (62 - yx 
44a;^) +2 (l + 6x + -a;M lna;-2a;ln(l -a;) - {A + x)\v? x + pfg{.x) 



17 3 
— + 7C(2)-31na;- - In^a; - 21na;ln(l - a;) - ln^(l - a;) +85i(a;) 

18 2 

+ Pfg{-x)S2{x)} , 
P^(^'J)(x) = C^Tn/ {-2 + 3a; + (-7 + 8x) Inx - 4 ln(l - x) + (1 - 2a;) In^ x + 2pGF(a:) 

-5 + 6C(2) + lnx-ln(l - x) - In^x - 21nxln(l - x) -ln^(l-x) 



+ SSi{x)\] + CATnf 

+ 2(1 + 4x) In^ X + pgf{x 



20 

9 K~~ 
178 

~9~ 



76 + 83x 



;i + 19x)lnx + 41n(l-x) 



14C(2) InxH ln(l - x) -In^x 



+ 81nxln(l-x)+21n2(l-x)-165i(x)] +2pgf{-x)S2{x)'] 



+ (^"-z)^ 1 -3 - -^Pgf{x} 



— h Inx + ln(l — x) 



P^i^ 



Ci ( 3 + 3C(3) j - CfTuj - -CATrif 



5{l-x) + 



C\[%- 2C(2) 
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20 „ ^ 



X 



CpTn 



23 
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7x H — x^ ) Inx + 2(1 + x) In x 



- -(1 + x) Inx - -Pgg{x) Inx 



9x 
CaTuj 

+ Cl|-^(25 + 109x)-2C(2) 



2( — + 5 
3x 



2 + X - x^ + 



X 



44 



X 



2/23 

+ 29 - 19x + 23x^ 

9 V X 



11 - 25x Inx -4(1 + x)ln X 



+ Pgg (a;) Inx 



22 



— 3 In X + 4 ln(l — x) 



+ 2pgg{-x)S2{x] 



(17) 



Here, C is Riemann's zeta function, with values C(2) = 7r^/6 and (^(3) ~ 1.202, and we 
have used the auxihary functions 

, , l + x2 
Pff{x) = , 

1 — X 
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, , l + (l-a;)^ 
Pfg{x) = , 

X 

Pgf{x) =x^ + {l-xf, 

[l-x(l-x)]^ 
Pgg[x) = 



x(l — x) 



and 



Si{x) = -Li2(l -x), 



S2{x) = -C(2) + -In^x - 21nxln(l + x) - 2Li2(-x), (19) 

where Li2(a;) = — Jq dt ln(l — t)/t is the dilogarithm. In the case of the colour gauge group 
SU(A^c), the Casimir operators of the fundamental and adjoint representations have the 
eigenvalues Cp = {N^ — 1) /{'^^Nc) and Ca = Nc, respectively. The trace normalization of 
the fundamental representation is T = 1/2. 
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TABLE CAPTIONS 

Table 1: Fit parameters for the charm- and bottom-quark FF's of sets S and M at LO and 
NLO. The corresponding starting scales are /Xq = 2mc = 3 GeV and /io = Sm^ = 10 GeV, 
respectively. All other FF's are taken to be zero at Hq = 2mc 

Table 2: x^ P^r degree of freedom pertaining to the LO and NLO fits of types S and M to 



the ALEPH |T^ and OPAL [^ samples of Z -^ cc and Z ^> hh and their combinations. 



The first bin of the OPAL data has been excluded 

Table 3: Branching fractions of charm and bottom quarks into D*"^ mesons at the respec- 
tive starting scales and at /^ = Mz evaluated from (1) with sets S and M at LO and NLO 



Table 4: Average momentum fractions of D*^ mesons produced through charm- and 
bottom-quark fragmentation at the respective starting scales and at /x = Mz evaluated 
from (|T0|) with sets S and M at LO and NLO 
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FIGURE CAPTIONS 

Figure 1: The cross sections of inclusive D*^ production in e^e~ annihilation evaluated 



with (a) set S and (b) set M at LO and NLO are compared with the ALEPH data |T5|. 
The three sets of curves and data correspond to the Z ^ cc and Z ^ bb samples as well 
as their combination 



Figure 2: As Fig. [l|, but for the OPAL data p 



Figure 3: The ARGUS |2g, HRS |T9[, TASSO [g, and OPAL 1|I§ data on D*± production 
in e~^e~ annihilation are compared with the LO and NLO predictions based on set S. For 
separation, the data have been rescaled by powers of 10. In the case of TASSO, the open 
triangles refer to the D^ — » K'^ti^ti^ti^ channel (TASSO 1) and the solid triangles to the 
D^ -^ K~TT~^ channel (TASSO 2). In the case of OPAL, we consider the dimensionless 
quantity {l/atot)d<j/dx and discard the data point at x = 0.05 

Figure 4: The LO and NLO predictions of inclusive D*^ photoproduction in ep collisions 
based on sets S and M are compared with the ZEUS data 0. We consider (a) the pr 
distribution da/dpT integrated over —1.5 < |/iab < 1 and 115 GeV < W < 280 GeV, (b) 
the j/iab distribution da/dyxs^^ integrated over 3 GeV < p^ < 12 GeV and 115 GeV <W < 
280 GeV, and (c) the W distribution da/dW integrated over 3 GeV < pr < 12 GeV and 
-1.5 < yiab < 1 

Figure 5: Influence of the bottom-quark FF's on the LO and NLO predictions of inclusive 
D*"^ photoproduction in ep collisions based on sets S and M. The analysis of Fig. ^a is 
repeated for Uf = 5 with and without the bottom-quark FF's, and the ratio of the two 
results is shown as a function of pt 

Figure 6: Influence of the charm-quark PDF of the photon on the NLO prediction of 
inclusive D*^ photoproduction in ep collisions based on set S. The analysis of Fig. ^ is 
repeated with this PDF switched off. For comparison, also the direct-photon contribution 
is shown 

Figure 7: Influence of the photon PDF's on the NLO prediction of inclusive D*^ photo- 
production in ep collisions based on set S. The analysis of Fig. lib is repeated using set 



GS HO [|r^ and the ACFGP ||T2| sets for massless and massive charm quarks 
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set 


flavour 


N 


a 


P 


Cc 


LO S 


c 


449 


8.01 


3.02 


- 


b 


335 


3.24 


6.47 


- 


NLO S 


c 


69.7 


8.78 


1.60 


- 


b 


260 


3.33 


5.83 


- 


LOM 


c 


0.202 


- 


- 


0.0856 


b 


235 


3.15 


6.06 


- 


NLO M 


c 


0.0677 


- 


- 


0.0204 


b 


214 


3.25 


5.65 


- 



Table 1 



set 


total 


sum 


\LEFI 

c 


I 
b 


sum 


OPAL 

c 


b 


LO S 


0.85 


0.84 


1.54 


0.78 


0.44 


1.05 


0.66 


NLO S 


0.90 


0.96 


1.61 


0.79 


0.52 


1.12 


0.55 


LOM 


1.05 


1.25 


2.43 


0.92 


0.37 


0.99 


0.68 


NLO M 


0.96 


1.21 


2.08 


0.82 


0.37 


0.94 


0.56 



Table 2 



set 


B,{2m,) 


BciMz) 


Bb{2mb) 


Bt{Mz) 


LOS 


0.220 


0.210 


0.242 


0.228 


NLO S 


0.217 


0.209 


0.236 


0.223 


LO M 


0.246 


0.231 


0.229 
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